Ground state of a confined Yukawa plasma including correlation eff'ects 
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The ground state of an externally confined one-component Yukawa plasma is derived analytically 
using the local density approximation (LDA). In particular, the radial density profile is computed. 
The results are compared with the recently obtained mean- field (MF) density profile il4|. While 
the MF results are more accurate for weak screening, LDA with correlations included yields the 
proper description for large screening. By comparison with first-principle simulations for three- 
dimensional spherical Yukawa crystals we demonstrate that both approximations complement each 
other. Together they accurately describe the density profile in the full range of screening parameters. 
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I. INTRODUCTION 

Interacting particles in confinement potentials are om- 
nipresent in nature and laboratory systems such as 
trapped ions, e.g. P, [1], dusty plasmas, e.g. [1, 0, Q 
or ultracold Bose and Fermi gases [1, An interest- 
ing aspect of particle traps is that it is easy to realize 
situations of strong correlations. The observed particle 
arrangements reach from gas-like, liquid-like to solid be- 
havior where the symmetry is influenced by the trap ge- 
ometry. Of particular recent interest have been spheri- 
cal traps in which plasma crystals consistin g of spherical 
shells (Yukawa balls) are formed, e.g. [sl. [ol llOl The 
particle distribution among the shells is by now well un- 
derstood [HEIil. 

In a recent study [l3| we have analyzed also the av- 
erage particle density in the trap and found that it is 
very sensitive to the binary interaction: it changes from 
a flat profile in case of long-range Coulomb interaction 
to a profile rapidly decaying away from the trap cen- 
ter in the case of a screened Yukawa potential. Using 
a non-local mean-field (MF) approximation the density 
profile could be computed analytically and was found to 
agree very well with first principle computer simulations 
for Yukawa crystals. However, when the screening is in- 
creased deviations in the trap center kept growing which 
were attributed to correlation effects missing in the mean- 
field model. 

The goal of this paper is to remove these discrepancies. 
For this we extend the analysis of Ref. [l3| by including 
correlation effects following an idea of Totsuji et al. [1^| 
applied to 2D systems. We apply the local density ap- 
proximation (LDA) using known results for the corre- 
lation energy of a homogeneous one-component Yukawa 
plasma. The results clearly confirm that correlation ef- 



fects are responsible for the strong density increase in the 
trap center. We find that LDA with correlations included 
agrees very well with simulations of Yukawa crystals in 
the limit of strong screening. On the other hand, for 
weak screening, the previous MF result turns out to be 
more accurate. Interestingly, for intermediate values of 
the screening parameter both methods are accurate, so a 
combination of both allows to quantitatively describe the 
density profile in the whole range of screening parame- 
ters. 

This paper is organized as follows. In Sec. II we in- 
troduce the LDA and use it first to compute the density 
profile in a mean-field approximation, which, of course, 
gives worse results than MF, but helps to understand 
LDA. Then, in Sec. Ill we improve the LDA model by 
including correlation effects. In Sec. IV the results for 
the density profile are compared to molecular dynamics 
simulations. A discussion is given in Sec. V. 



II. GROUND STATE OF A CONFINED 
PLASMA WITHIN LDA 

We consider N identical particles with mass m and 
charge Q confined by an external potential $ and in- 
teracting with the isotropic Yukawa-type pair potential, 
^(^) — {Q"^ /t) exp(— Kr). To derive the properties of in- 
terest we start with the expression of the ground-state 
energy, which is given by 



E[n] 



(1) 



with the energy density u(r) = M™^(r) -I- u'^°'^{r), where 
the mean-field energy density is 



u'"'(r) =n(r)<^ $(r) -|- 



A^- 1 
2N 



J d\2n{r2)V{\r-r2\)y 
(2) 
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(Sec. Ill) by means of the local density approximation. 
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Before, we introduce this approximation and obtain first 
the LDA results in mean-field approximation (LDA-MF). 
These results will not be as accurate as the MF results, 
due to the applied approximation, but LDA-MF helps to 
familiarize oneself with LDA and its characteristics. 

The local density approximation is based upon the idea 
of replacing the non-local terms within the energy den- 
sity at point r by local expressions using the known en- 
ergy density of the homogeneous system with its density 
no equal to the local density n(r) of the true inhomoge- 
neous system in question. Therefore, to derive LDA-MF 
we need to substitute for the non-local second term in 
([2]), i.e. for the density of interaction energy, the corre- 
sponding expression of the infinite homogeneous system, 
which is given by [details are given in the Appendix] 
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uo{k) = no 
N 



-^Q d\.no- 



1-2 I 



(3) 



TV ^ ""^2' 



and, as a second step, replace the homogeneous density 
no by the local density n(r). Thus we obtain the LDA- 
MF ground-state energy 



mf r 

ldaI 



d^r-u(r) 



(4) 



with the energy density 



u{v) =n(r)|$(r)-f ^^^Q2^(r)J 
The variation of the energy 

with respect to the density 7i(r) (for details see Ref. 
yields an explicit expression for the density profile in an 
arbitrary confinement potential 



(5) 



(6) 



n(r) = — ^ — TT U - $(r) 

^ ' 47r(7V- 1)02 ^ 



(7) 



which holds at any point where the density is nonzero. 
Due to ^ this density is normalized by 



(frn(r) = N. 



(8) 



The case of an isotropic confinement $(r) = $(r), 
which is of particular interest, leads to an isotropic den- 
sity distribution n(r) = n{r) = n{r)Q{B. — r) the outer 
radius R of which is being fixed by the normalization 
condition ([5]) which now becomes dr r^n{r) = N/At:. 
In this isotropic case the yet unknown Lagrange multi- 
plier /i can be obtained by taking the variation also with 
respect to R , which yields 



ii = ^{R). 



(9) 
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FIG. 1: Radial density profile for a parabolic confinement 
potential <[>(r) = (a/2)r^ and a constant screening param- 
eter K,dc = 1 and three different particle numbers A'^ = 
100, 700, 2000. Also the result for K,da = 0.4, N = 2000 is 
shown by the dashed line. For comparison, the non-local MF 
results for Kdc = 0.4, 1.0, N — 2000 are given by the dotted 
lines. 



Compared to the MF result which was given in [1 
N 



n"^'(r) 



47r(A^ - l)g2 



(10) 
(11) 



the LDA-MF density ^ differs in two points. On the 
one hand the Laplacian of the potential A$(r) is missing 
and, on the other hand, the expression of the chemical 
potential /i is simpler than /i™^. That is based upon the 
fact, that the missing expressions consist of derivatives 
and thus contain informations about contiguous values 
of the potential, what is suppressed within LDA-MF and 
generally within LDA. 



A. Parabolic confinement potential 

For the case of a parabolic external potential $(r) 
{a/2)r'^ the density following from Eqs. ([7]) and ([9|) is 



n{r) 



aN 



47r(iV - 1)Q2 



Q{R-r). (12) 



The dimensionless combination k_R, which contains the 
limiting outer radius, can be obtained from the normal- 
ization ([8]) and is given by 



kR 



5/l5(iV- 1)02^3 



a 



J 15 



(13) 



(K4)3(iV-l). 
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FIG. 2: (Color online) Radial MF density profile (solid lines) 
compared to LDA-MF (dashed lines) for three diflterent den- 
sity parameters x — (ndc)'^ {N — 1). The abscissa is normalized 
with the MF radius i?™' , while the ordinate is normalized with 
the corresponding MF density n™'(0) at the trap center. 



Here, we introduced the length scale c?c — (SQ^/a)^^'^, 
which is the stable distance of two charged particles in 
the absence of screening 11] and which will be used below 
as the proper unit for lengths and screening parameters. 
As unit for densities we use the average density of a large 
Coulomb system, which is given by Uc = {3a) / (AttQ'^) . 

The results of (fT2|) are shown in Fig. [1] for three par- 
ticle numbers from N = 100 to = 2000. One clearly 
sees the parabolic decrease of the density away from the 
trap center till it terminates in zero. The curvature of 
the density does not change by increasing the particle 
number - just the density increases continuously at every 
space point and, at the same time, extends to higher val- 
ues of the limiting radius R. However, the curvature of 
the density profile changes dramatically, when the plasma 
screening is increased at constant N. 

Thus, in the case of an isotropic parabolic potential, 
the LDA density profile bears qualitative resemblance to 
the density profile in the non-local mean-field approxima- 
tion, although in the case of other confinement potentials 
the deviations of the LDA-MF from the MF approxi- 
mation are stronger, cf. Eqs. ([7]) and pU)) . However, 
quantitatively in two points MF differs from LDA-MF 
for parabolic confinement as well, as can also be seen in 
Fig-Si First, the density in this local density approxima- 
tion does not show a discontinuity at r = i?, in contrast 
to the MF resuh, Eqs. (dHl E]). This is due to the ne- 
glect of edge effects in this derivation of the LDA result. 
Secondly, LDA-MF yields too small values of the density. 
This error is reduced (cf. Fig. (2) with increasing values 
of the density parameter x = (ndc)^ {N — 1), cf. Ref. p^ . 
which regardless of the factor N/{N—l) solely determines 
the density profile. The reason for this improved behav- 
ior with increasing x is due to the fact that an increase 
of K contracts the effective area of integration within ^ 
as well as within ([3]). The contraction finally is in favor 



of the accuracy of LDA-MF, because the decreased in- 
tegration volume contains a more homogeneous density. 
Also an increase of the particle number N, what flattens 
the density profile, will similarly improve LDA-MF. 

Because the validity of the mean-field model depends 
on the value of the screening parameter ndc, there are 
the following two cases. In the first case, for small val- 
ues of the screening parameter, the MF approximation 
provides a good description of the density profile, but 
LDA-MF underrates this profile and so does not give a 
good description on its own. (That applies also if finite- 
size effects are included, cf. Fig. [3l) In the second case, 
for large values of the screening parameter, the LDA-MF 
approaches the MF approximation, however, there, the 
latter does not describe the density profile correctly due 
to the neglect of the now relevant correlation contribu- 
tions p^ . Thus, the local density approximation of the 
mean-field energy alone does not give a good description 
of the density profile. 

However, it gives a straightforward way to include the 
missing correlation contributions in the energy density 
by usage of the result of the homogeneous system, see 
Sec. imi 



B. Improvement of LDA by inclusion of finite-size 
effects 

As can be seen from Fig. [5] and from Eq. ^ the den- 
sity profile obtained by LDA-MF breaks down in the 
Coulomb case - the density cannot be normalized any- 
more, which is the same as in the two-dimensional case 
[l5| . But the application of a local density approxima- 
tion cannot be the reason for this, because the method of 
LDA is based upon the usage of results from the homoge- 
neous system, and the Coulomb system is homogeneous 
with no = [N/{N - l)]nc. 

In fact, the cause of the breakdown is the use of results 
from the infinite homogeneous system neglecting finite- 
size effects. This failure can be avoided by replacing ([3]) 
by the corresponding expression of the finite homoge- 
neous system. In the appendix such an expression is de- 
rived for isotropic confinement. As a result the finite-size 
effects lead to a corrected density profile 



n(r) 



47r(iV- 1)Q2 



^ - <I>(r) 



1 — e~'*-^(l + kR) sinh(Kr)/(Kr) 



e(i?-r), 



(14) 



instead of Eq. (|12p . what indeed yields the constant 
MF solution also for LDA-MF. As another example, in 
Fig. [3] the density profiles with [LDA-MF (fs corrected)] 
and without these finite-size contributions are shown for 
N — 1000, ndc = 0.3. One clearly sees, that in the 
case of finite-size correction the density profile shows a 
discontinuity at the boundary and, due to that, it yields 
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FIG. 3: Radial density profiles of a spherical plasma of 
iV = 1000 and Kdc = 0.3 calculated by LDA-MF with (fs 
corrected) and without finite-size effects included. For com- 
parison the exact MF result is also given by the dashed line. 
The difference between the finite-size correction and the par- 
tial finite-size correction is described in the text. 



increased values of the density. However, the density pro- 
file including edge effects is not monotonically decreasing 
away from the trap center but has a density increasing 
part in the outer range, what is not correct. This is due 
to the spatial dependence of the denominator of Eq. 

By contrast a more accurate monotonically decreasing 
density profile can also be obtained by taking the finite- 
size effects only partly into account [LDA-MF (partly fs 
corrected)], as derived in the appendix. The final result 
is given by 



n(r) 



fi ~ <i>(r) 



47r(A^ - 1)Q2 1 _ e-«-R(l + ^R) 



-e(i?-r), (15) 



which now misses the spatial dependence in the denomi- 
nator. The corresponding result is also given in Fig. [3) 

Consequently, for Yukawa systems like those analyzed 
here an improvement of LDA is possible by including 
edge effects. However, for small values of the screening 
parameter even the improved local density approxima- 
tion does not approach the degree of accuracy obtained 
by the non-local mean-field model, cf. Fig. [31 On the 
other hand, for increased screening the finite-size effects 
do not alter the density profile significantly. 

Therefore, below we continue to use Eq. ([3]) of the in- 
finite homogeneous system, what will not interfere the 
following results. 



III. INCLUSION OF CORRELATION 
CONTRIBUTIONS 

The energy expression i?LDA (HE]) contains only the 
energy density of the confinement and of the mean-field 
interaction. To include the contribution of the particle 
correlations we can make use of the result for the density 



of correlation energy of the homogeneous system which 
is given by Eq. (3) of Ref. 



Ucorr(?^o, «:) = - 1.444Q^n^ 

X exp|^— 0.375KnQ 



7.4-10" 



(16) 



where tt-q is the corresponding density of the homoge- 
neous system. By replacing this density with the local 
density n(r) of the inhomogeneous system one obtains 
the correlation contribution of the energy density within 
LDA. Thus we derive the complete ground-state energy 
in local density approximation 



(fir u{r) 



(17) 



with energy density 

u(r) =n(r)$(r) + ^^^tt^Q^ n(r)^^ - lAAAQ^nM^ 

X exp(^-0.375Kn(r)"3 + 7.4 • 10~'^(Kn(r)"3)'') . 

(18) 

As before, variation of the energy ([T7)) at constant parti- 
cle number, cf. Eq. yields the ground state density 
profile, but now with correlation effects included. In this 
case the strong non-linear character of the energy density 
does not allow for an explicit solution. Just an implicit 
solution is possible and is given by the following equation 
for z'^(r) = K~^n{r), which can be regarded as the local 
plasma parameter of the system, 

=^^^(r) + - (ciz(r) + c,- c,zir)-') 

X exp(-0.3752(r)"i + 7.4 • 10"^z(r)" 

The constants Ci are given by 



ci 

C3 



0.153 
0.0144 
1.134- 10-^ 



(19) 

(20a) 
(20b) 
(20c) 



The solution of Eq. (1191) '^^^ be obtained numerically. 
For the case of a parabolic external potential $(r) = 
{a/2)r'^ results are given in Fig. 21 There, the density 
profile of a plasma of iV = 2000 particles within LDA 
is shown for three different screening parameters: ndc — 
0.5, Kdc = 2.0 and ndc — 3.0. For comparison the LDA- 
MF density profile is shown too. 

It can be seen that for a small screening parameter (see 
line ndc = 0.5) both density profiles are nearly identical. 
But with increasing screening, i.e. for smaller values of 
the local plasma parameter z'^, the correlation contribu- 
tions within LDA alter the curvature of the profile, which 
rises more steeply towards the center. So the particle 
correlations tend to increase the central density of the 
plasma, which also can be seen in Fig. O in comparison 
with the mean-field approximation. 
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FIG. 4: Radial density profile of a confined spfierical plasma 
of N = 2000 particles calculated with LDA including correla- 
tion contributions (solid lines) compared to LDA-MF (dashed 
lines) for three different screening parameters. 



IV. COMPARISON WITH SIMULATION 
RESULTS FOR FINITE YUKAWA CRYSTALS 



We performed molecular dynamics simulations of the 
ground state of a large number of Coulomb balls for the 
purpose of comparison with the average density of spher- 
ical Yukawa crystals, see refs. [13, Ell for details. In order 
to obtain a smooth average radial density profile the av- 
eraging process was accomplished by substituting each 
particle by a small but finite sphere. Numerical results 
of this comparison with a Coulomb ball of iV = 1000 
particles are included in Fig. [5] for four different screen- 
ing parameters. The symbols denote the average particle 
density in the vicinity of the corresponding shell, while 
the lines represent the MF (solid) and the LDA density 
(dashed). 
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FIG. 5: Radial density profiles of a three-dimensional plasma 
of A'^ = 1000 particles calculated with the exact mean-field 
model (solid lines) and with LDA including correlation con- 
tributions (dashed lines) for four different screening parame- 
ters: Kdc = 1, ndc = 2, ndc = 3 and ndc — 5. Averaged shell 
densities of molecular dynamics results of a plasma crystal for 
the same parameters are shown by the symbols. For details 
see discussion in Sec. IIVI 



SUMMARY AND DISCUSSION 



For small values of the screening parameter ndc < 2 the 
simulation results are very well reproduced by the ana- 
lytical density profile of the non-local mean-field model 
(MF), whereas the local density approximation under- 
rates the results (lower lines in Fig. [5] (a)). On the other 
hand, for larger values of the screening parameter kcIc > 2 
the simulation results are reproduced by LDA, whereas 
MF underestimates these results in the center. This un- 
derestimation is accompanied by a wrong prediction of 
the profile curvature (Fig. O (b)). For intermediate val- 
ues of the screening parameter ndc ~ 2 both methods 
are very close to the averaged simulation results (upper 
lines in Fig. \5\ (a)). We have verified this behavior also 
for other Coulomb balls. Another representative example 
is shown in Fig. [S]for a Coulomb ball with N — 10000. 
There, the same behavior as in Fig. [5] is seen. 



A theoretical analysis of the ground state density pro- 
file of a spatially confined one-component plasma within 
local density approximation was presented. We derived 
a closed equation, Eq. for the density profile in- 

cluding correlations effects for arbitrary confinement po- 
tentials with any symmetry. In contrast to the result 
without particle correlations the density profile shows an 
increased central density with increasing screening pa- 
rameter. The validity of LDA is however limited to not 
too small values of the screening parameter, Kdc > 2. 

Comparisons with first-principle simulation results of 
strongly correlated Coulomb clusters with varying screen- 
ing parameter showed that LDA allows to remove the 
problem of the MF approximation observed in Ref. [l3| 
which arises with increasing screening parameter. There- 
fore, the mean-field model together with the presented 
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FIG. 6: Radial density profiles of a three-dimensional plasma 
with A'^ — 10000 and two different screening parameters 
{hide = 1, Kdc = 3). The solid (dashed) lines show MF (LDA) 
calculation results. Symbols denote molecular dynamics re- 
sults of a plasma crystal for the same parameters where the 
average density at the positions of the shells is shown. 



local density approximation complement one another in 
the description of strongly correlated spatially confined 
one-component plasmas. 



This formula has no spatial dependence due to the in- 
finite integration volume and it diverges in the limit of 
Coulomb interaction (k 0) leading to a breakdown of 
the approximation. 

By contrast, the density of interaction energy of the 
corresponding finite homogeneous system (a sphere with 
center r2 = and radius R) is given by 



uo(K,r) 



N - 1 



no 



2N 



(fir 2 no- 



,-K|r-r2| 

Ir-ral 



g — e 



2N KT 



^> n. 



dr2 r2 sinh(Kr2) 



N 



+ sinh(Kr) / dr2 '^^^ 

J r 

l-e-««(l + .i?)^^^ 



(A2) 
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APPENDIX A: LOCAL DENSITY 
APPROXIMATION USING A FINITE 
REFERENCE SYSTEM 



including a finite-size contribution, which prevents the 
problem of divergence at k — > 0. As already mentioned in 
Sec. Ill BI the resulting density profiles show the incorrect 
behavior of a non-monotonic density profile, cf. Fig. [31 

An improved correction which only partly takes edge 
effects into account can yet be obtained by removing the 
explicit spatial dependence within the density of interac- 
tion energy 



The investigation of an inhomogeneous system within 
LDA uses known results from the corresponding homo- 
geneous system. There, the infinite homogeneous system 
is often used as reference system what entails that finite- 
size effects are neglected. To take these into account the 
finite homogeneous system has to be used as reference. 
In the present derivation such a modification is made for 
an isotropic confinement and leads to a change of the 
expression for the density of interaction energy, Eq. ([3]), 



rR 



,7V- 1 



n, 



N 

N - 1 
N 



Q 2t: dr2 r2e 



,2tt 



Q^nl—{l-e-''"{l + KR)) 



(A3) 



wo(k) = no 



N - 1 
2A^ 



Q2 / (fr2no- 



2N 



-K|r-r2| 



r - r2| 

— Kr2 



r2 



(Al) 



This expression also has no divergent limit for k 0, 
and, at the same time, yields monotonically decreasing 
density profiles as can also be seen in Fig. [3] 
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